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> . ABSTRACT 

J^j ' We present broadband multi-wavelength observations of GRB 080310 at redshift z = 

r> \ 2.43. This burst was bright and long-lived, and unusual in having extensive optical 

^ . and near IR follow-up during the prompt phase. Using these data we attempt to 

simultaneously model the 7-ray, X-ray, optical and IR emission using a series of prompt 
pulses and an afterglow component. Initial attempts to extrapolate the high energy 
model directly to lower energies for each pulse reveal that a spectral break is required 
between the optical regime and 0.3 keV to avoid over predicting the optical flux. We 
demonstrate that afterglow emission alone is insufficient to describe all morphology 
seen in the optical and IR data. Allowing the prompt component to dominate the 
early-time optical and IR and permitting each pulse to have an independent low energy 
spectral indices we produce an alternative scenario which better describes the optical 
light curve. This, however, does not describe the spectral shape of GRB 080310 at 
early times. The fit statistics for the prompt and afterglow dominated models are 
nearly identical making it difficult to favour either. However one enduring result is 
that both models require a low energy spectral index consistent with self absorption 
for at least some of the pulses identified in the high energy emission model. 

Key words: gamma-rays: bursts. 
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1 INTRODUCTION 

Over the last few years a combination of fast-response 
ground-based telescopes triggered by the availability of 
rapid, accurate localisations have started to provide the data 
required to answer the question of what is causing the early, 
bright X-ray and optical emission from gamma-ray bursts 
(GRBs). The most acc urate prompt X-ray locations come 
from the Swift satellite (|Gehrels et alj|2004 ). These are sup- 
plemented by either on-board or ground detections of the 
ultraviolet (UV), optical or infra-red (IR) counterpart. 

In the popular relativistic fireball model for GRBs, the 
early, usually highly variable emission is understood to be 
due to internal shocks (I Sari fc Piranlll997l ) or magnetic dissi- 
pation within the jet, and the so-called external emission is 
produced by the interaction of the jet with the surrounding 
medium. The l atter emission is usually described using the 
fireball model (|Rees fc Meszaro^ [l992). which has success- 
fully been applied to describe the behaviour of GRBs half 
a day or so after the trigger, but has difficulties explaining 
the complex behaviour seen in the first few hours, a period 
now routinely accessed by Swift and other rapid-response fa- 
cilities. Ideally multi-wavelength observations should be ob- 
tained while the burst is happening so as to try to disentan- 
gle the relative contribution from the internal and external 
com ponents. 

lEvans et alj (|2009l ) present a uniformly analysed com- 
prehensive sample of 317 Swift GRBs spanning from De- 
cember 2004 to July 2008, in which their morphologies 
are compared to th e proposed canonic al X-ray light curve 
dNousek et all I2006L IZhang et all 120061 & iPanaitescu et all 
120061 ). Such canonical light curves consider the X-ray emis- 
sion to consist of a series of power-laws, where one impor- 
tant phase is the rapid decay phase which has been ex- 
plained as being the smooth continuation of th e prom pt 
emission JO'Brien et ail |200d . iTagliaferri et al.l 120051 fc 
iBarthelmv et al.l I2005T ) . From the sample of lEvans et al.l 
(2009) it is clear that the X-ray light curves of GRBs vary 
from burst to burst. Some show stro ng flaring, one example 
being GRB 061121 (|Page et al.ll2007h , however, some bursts 
show remarkably s imple and smooth decay (GRB 061007; 
Schad v et al.l 20 07). Similar findings are also reported in 



Racusin et alj (|2009T ). Rapid behaviour, such as flaring, at 



high energies is often attributed to central-engine behaviour 
l|Margutti et alj l20ld ) but how this relates to the optical 
emission remains somewhat of a mystery. The available 
datasets reveal a confusing picture. In some cases the early 
optical data seem to tra ce the X-ray and 7-ra y light curves, 
such as GRB 041219A (|Vestrand et al.ll2005l fc iBlake et all 
120051 ) , suggesting that optical flaring may be of internal ori- 
gin. In other GRBs the optical behaviour seems entirely 
unrelated to the t he high-energy emission (GRB 990123; 
lAkerlof et a0l999lfc G RB 060607A; lNvsewander et"aTll2009l 



& Molinari et al 



2003), and instead seems to follow the be- 



haviour of the external afterglow. 

To make progress requires continued efforts to observe 
GRBs over as wide a wavelength range as possible from 
as early as possible. This is only really viable for bright, 
long-lasting GRBs which are well-placed for rapid follow- 
up. Here we present prompt, multi-wavelength data from the 

f E-mail: oml2@star.le.ac.uk (OML) 



GRB 080310, which begin at a time soon after the trigger 
that is not often accessed by ground-based facilities. Follow- 
ing the trigger, GRB 080310 was det ected on-board Swif t 



by both the X-Ray Telescope (XRT) (Burrows ct al 



2005) 



2005) 



and UV/Optical Telescope (UVOT) (|Roming et al 
and also observed in the optical and IR by several ground- 
based telescopes. These rare data present us with an op- 
portunity to discriminate between whether the early time 
lower-energy light curve of a GRB is driven by internal or 
external emission, during which time the high-energy emis- 
sion is presumed to be totally internally dominated. 

In the following section we discuss the observations, 
then in §3 we present attempt s to fit both an internal shock 
model dGenet fc Granotl [2009) and an afterglow component 
l|Willingale et all I2007T ) to the X-ray and 7-ray data, and 
describe the necessary modification required to simultane- 
ously fit the early optical emission in two scenarios, where 
either the prompt or afterglow components dominate this 
early flux. In §4 we briefly discuss the relative merits of us- 
ing each component, before finally in §5 we conclude which 
of the two alternatives is a better fit to the observed data 
and consider the implications of the model on the physics 
governing the emission from GRB 080310. 



2 OBSERVATIONS 

On 2008 March 10 the Swift BAT triggered and lo- 
cated GRB 080310 (trigger numbe r 305288) on board at 
08:37:58 UT jCummings et alj|2008h . Swift slewed immedi- 
ately which enabled the narrow field instruments to begin 
observing the burst 89 s after the trigger. The burst was 
detected by the XRT and UVOT (white filter), with the 
latter providing the best Swift position of RA(J2000) = 
14 fc 40 m 13 s g0) Dec(J2000) = -00°10'29'.'6 with a la error 
radius of 0'.'6. Figure [1] shows a UVOT f-band image from 
the early time data. The Swift light curves obtained in mul- 
tiple bands from each of the three on board instruments 
are presented in Figure [2] and most of the available datasets 
from an extensive number of facilities are show n in Figure 
El Data from the PAIRITEL jPerlev et alj|2008l ) and KAIT 
instruments are shown separately in Figure [4] 

In addition to Swift observations, and those instruments 
already mentioned, GRB 080310 was also observed on the 
ground with numerous op tical and NIR facilities, inc luding 
REM |Covino et alj l2008h . VLT (|Covino et all l200St) and 
the Faulkes Telescope North. These observations are shown 
alongside the BAT and XRT light curves in Figure [3] 

The Kast dual spectrometer at the Lick Observatory, 
California, obtained t he first redshift estim ation for this 
burst of z = 2.4266 (|Prochaska et al.1 120081 ) using strong 
absorption features from Silicon, Carbon and Aluminium. 
This was later corroborated by the VLT-UVES instrument 



Vreeswiik et al. 


2008) 


Prochaska et al. 


2008) 



The following subsections describe the observations in 
more detail. The Swift data analysis was performed using re- 
lease 2.8 of the Swift software tools. Parameter uncertainties 
are estimated at the 90% confidence level. We note that the 
optical datasets have been reduced using different methods, 
and fully investigated the effects of cross calibration errors, 
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Figure 1. UVOT D-band image from Tp +79 s to T + -1539 s. The 
enhanced Swift XRT position is shown ( Osbor ne et al . 2008), with 
an error circle of radius l'/4, corresponding to the 90% confidence 
limit, containing the optical afterglow candidate at an average 
■u-band magnitude of 17.49. 



ensuring that our later analysis remained insensitive to their 
effects. 



2.1 BAT 

The BAT data were processed using the standard bat- 
grbproduct script. The top four panels of Figure[2]show the 
BAT light curves displayed in the standard energy bands of 
15 - 25keV, 25 - 50keV, 50 - 100 keV and 100 - 150 keV, 
plotted with respect to the BAT trigger time (To)- The bin- 
ning is such that each bin satisfies a minimum signal-to-noise 
ratio of 5 and a minimum time bin size of 1 s. 

The 7-ray light curve shows many peaks with the first at 
To - 60s. The brightest peak extends from To — 12 s to To+7s. 
This is followed by a period of no detectable emission before 
a weaker, broad series of pea ks is seen from To + 180 s to 
T + 360 s (|Tueller et alJ lioOS - ). The latter peak is consistent 
with the first strong flare seen in the XRT (see below). The 
BAT emission is strongest in the lower energy bands, below 
~ 100 keV. The T 90 is estimated to be 365 ± 20 s (where the 
error includes systematics). 

The total spectrum from T - 71.76 s to T + 318.75 s is 
well fit by a power-law of photon index 2.32 ± 0.16, with a 
total fluence of 2.3±0.2 x 10" 6 erg cm" 2 over the 15- 150 keV 
band. The fluence ratio S(25 - 50keV)/S(50 - 100 keV) is 
1.27 ± 0.17 which puts GRB 080310 on the border of the X- 
ray-rich gamma-ra y bursts and X-ray flash (XRF) according 
to the definition of lSakamoto et aL I (|2008l ). 



2.2 XRT 

The Swift XRT started observations of GRB 080310 89 s 
after the trigger, with windowed timing data ranging from 



0.04 
0.02 

o.oo 

0.04 

0.02 

0.00 
0.04 

0.02 

0.00 
I 0.04 



S 0.00 
80.0 

40,0 

0.0 
80.0 

40.0 

_n 
f °' 3 

J 0.4 
" 0.0 



BAT 100-150 ke\i 



BAT 50-100 keV 



BAT 25-50 keV 



BAT 15-25 keV 





XRT 0.3-1.5 keV 




-100 100 200 300 400 500 600 700 800 
time since trigger (s) 



Figure 2. Swift prompt light curves from GRB 080310 in dif- 
ferent energy bands. From top to bottom, BAT 100 — 150 keV, 
BAT 50 - 100 keV, BAT 25 - 50keV, BAT 15 - 25keV, XRT 
1.5 - 10.0 keV, XRT 0.3 - 1.5 keV and UVOT v and white bands. 
Strictly, the count rates for BAT are per fully illuminated detec- 
tor. 



To + 95 s to Tp + 799 s and ph oton counting mode data after 
this l|Beardmore et aDl200&t ). As can be seen in Figure [2] 
the soft band XRT light curve initial fades until To + 130 s, 
before the burst rebrightens. This flaring activity continues 
until To + 420 s before the count rate drops again briefly. A 
further flaring event is seen between To + 500 s and To + 620 s, 
which is approximately half as bright as the first. There is 
significant spectral evolution during the flaring events, where 
the S(1.5- 10 keV)/S(0.3- 1.5 keV) ratio shows a hardening 
in the spectrum at T + 135 s to T + 200 s and T + 500 s to 
To + 565 s before, in both instances, there is a softening as 
the flaring behaviour declines. 



2.3 UVOT 

UVOT observations began with a 100 second finding char t 
exposure taken at To + 99 s l|Hoversten fc CummingsTl2008h . 
This finding chart exposure was the sum of 10 individual 10 
second exposures, with GRB 080310 being detected in all 
but the first exposure. At the time of GRB 080310 UVOT 
observations during the first orbit were taken in the event 
mode, which allows for higher timing resolution, while dur- 
ing subsequent orbits observations were taken in imaging 
mode. UVOT photometry was done using the publicly avail- 
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Figure 3. 7-ray to near IR light curves for GRB 080310, showing the energy bands used to observe the GRB and the instruments used 
to take the measurements. Data have been scaled from their actual values for viewing purposes. 



able FTOOLS data reduction suite, and is in th e UVOT 
photometric system described in lPoole et al.l(|200ct ). During 
the first 1,000 seconds the light curve is complex as shown 
in Figure [2] after which the the burst can be seen to fade in 
all the observed optical bands (Figure [3}. 

There are available data for the UVOT white, v, b, u 
and uvwl. As this burst has a measured redshift of z = 2.43 
it is important to consider absorption from the intervening 
medium through which the emission must travel. Correc- 
tion factors were calculat ed for Lyman a bsorp tion using the 
GRBz code described i n ICurran et al. (2008). GRBz uses 
the model presented by iMadaul (jl995h to calculate the ab- 
sorption from neutral Hydrogen in the intergalactic medium. 
Having found the correction factor for all the optical and 
near infrared bands presented in this work, we noted that it 
was only necessary to correct the data in the b, u and uvwl 
bands. Whilst the correction to both the u and b bands was 
not large, we found that the uvwl data required a signifi- 
cant correction, within which there was a large uncertainty. 
Given that there were only two data points from UVOT 
in this band, we decided to remove them from any further 
analysis. 



2.4 Faulkes Telescope North 

Observations with the Faulkes Telescope North (FTN) 
started at 2008-03-10 09:31:07.3 (UT), 3.188 ks after the 
trigger. Data were reduced in a standard fashion using iraf 
l|Todvlll986l ). Cali bration was perform ed using the SDSS 
data for the region |Adelman-McCarthv et al.|[2007l ). For the 
I and R filters, the SDSS photometry was converted to the 



Johnson-Cousins systemQ Photometry was then performed 
using an aperture matched to the average seeing of the 
(combined) frames. For the conversion from magnitude to 
flux, the data were first corrected for Galactic extinction us- 
ing th e COBE-DIRBE extinction maps from ISchlegel et al.l 
(1998 ), and then con verted using flu x zero points from 
iFukugita et al.l (|l995l ) for optical and iTokunaga &: Vaccal 
( 2005 ) for infrared. AB magnitudes were converted following 
lOke fc Gunnl (|l983l ). 



2.5 ROTSE 

ROTSE-IIIb, located at McDonald Observatory, Texas, re- 
sponded to GRB 080310 and began imaging 5.7 seconds af- 
ter the GCN notice time (|Yuan et al.ll200£ ). Observations 
were carried out under fluctuating weather conditions. The 
optical transient (OT) was detected between 25 minutes and 
3.5 hours after the trigger. To improve detection signal to 
noise ratio, sets of 4 to 11 images are co-added and expo- 
sures badly affected by weather are excluded. The OT is 
slightly blended with two nearby stars in the ROTSE im- 
ages. We therefore subtract the scaled point spread func- 
tions (PSFs) of these two nearby stars and then apply PSF- 
matching photometry on the OT using our custom RPHOT 
package (|Quimbv et al.|[2lll)(i). The analysis is further com- 
plicated by large seeing variation, particularly towards the 
end of the observation. The structure seen in the light curve 



1 See Lupton (2005) at 

http://www.sdss.org/dr6/algorithms/sdssUBVRITransform.html 



The early time emission of GRB 080310 5 



during this time is likely not significant. The ROTSE-III un- 
filtered magnitudes are calibrated to SDS S r using standar d 
stars in the pre-burst SDSS observations l|Cool et al.ll2008h . 



2.6 REM 

Early time optical and NIR data were collected using the 60- 
cm r obotic telescope REM (|Zerbi et alj|2003 ; ICovino et al.l 
120041 ') located at the European Southern Observatory (ESO) 
La Silla observatory (Chile). The telescope simultaneously 
feeds, by means of a dichroi c, the two focal instruments: 
REMIR (jConconi et al J |2004 > a NIR camera, operating in 
the range 1.0 to 2.3 /im (z' , J H and K) and R OSS (REM 
Optical Slitless Spectrograph; iTosti et alj|2004 ) an optical 
imager with spectroscopic (slitless) and photometric capa- 
bilities (V, R, I). Both cameras have a field of view of 10 x 
10 arcmin 2 . 

REM reacted automatically after receiving the Swift 
alert for GRB 080310, and began observing about 1 50 sec- 
onds after the GRB trigger time (|Covino et al.|[2008l ). 

Optical and NIR data were reduced following standard 
procedures. In particular, each single NIR observation with 
REMIR was performed with a dithering sequence of five 
images shifted by a few arcseconds. These images are au- 
tomatically elabor ated using the proprietary routine AQuA 
l|Testa et al.ll2004T ). The script aligns the images and co-adds 
all the frames to obtain one average image for each sequence. 
Astrometry was performed using USNO-BlQ and 2MAS93 
catalogue reference stars. 

Photometry was derived by a comb ination of the SEx- 
tractor package l|Bertin fc Arn outs 1996]) and the photomet- 
ric tools provided by the GAli^j package. The photometric 
calibration for the NIR was accomplished by applying aver- 
age magnitude shifts computed using 2MASS isolated and 
non-saturated stars. The optical data were calibrated using 
instrumental zero points and checked with observati ons of 
standard stars in the field provided bv lHendenl (|2008l ). 



the tools provided by the GAIA package and with PS F pho- 
tome try with the ESO-midas0DAOPHOT context l|Stetsonl 
Il987l ). 

The photometric calibration for the NIR was accom- 
plished by applying average magnitude shifts computed us- 
ing 2MASS isolated and non-saturated stars. The optical 
data were calibrated using instrumental zero points and 
checked with observatio ns of standard stars in the field pro- 
vided by iHendenl (|2008T ). Linear polarimetry position angle 
was corrected by means of observations of polarimetric stan- 
dard stars in the NGC 2024 region. 



2.8 WHT 

Late imaging was obtained with the 4.2m William Herschel 
Telescope (WHT), at Roque de los Muchachos Observatory 
(La Palma, Spain) using the Long-slit Intermediate Resolu- 
tion Infrared Spectrograph (LIRIS) in its imaging mode. Ob- 
servations consisted of 36 x 25 second exposures in i/-band, 
obtained on the 14th March 2008 from 04:36:28 to 04:53:41 
UT. The data were reduced following standard procedures 
in iraf. For the photometric calibration, we used stars from 
the 2MASS catalogue as reference. 



2.9 PAIRITEL 



PAIRITEL (|Bloom et al.ll2006l ) responded to GRB 080310 
and began taking data at 09:04:58 (UT) in the J, H, 
and K filters simultaneously (|Perlev et al.ll2008l ). The after- 
glow (IChornock et al. 20081 ) was well-detected in all three 
filters. Perlev et al. ( 20081 ) also report on an SED con- 



structed using data from PA IRITEL, KAIT and UVOT 
l|Hoversten fc Cummingsl 120081 ). allowing a joint fit to be 
made and the estimation of a small amount (Av ~ 0.10 
±0.05) of SMC-like host-galaxy extinction. 



2.7 VLT 

VLT FORS1 V and R observations for GRB 080310 were 
automatically activated with the RRM mod^fl allowing the 
telescope to react promptly to any alert. The field was ac- 
quired and the observations began less than seven minutes 
after the GRB trigger. Later VLT observations were ob- 
tained with ISAAC at about one day after the burst with 
the J, H and K filters. In addition linear polarimetry ob- 
servations were carried out with FORS1 with the V filter at 
approximately one, two and three days after the trigger. 

Optical and NIR data were reduced following stan- 
dard procedures_with the tools of the ESO-eclipse pack- 
age (|Devillardlll9"97l ). Polarimetric data were reduced again 
follow i ng standard procedures as discussed in ICovino et al.l 
l|l999l.l2002| . |2003l ). Photometry was performed by means of 



2 http:/ /tdc-www. harvard.edu/catalogs/ubl. html 

3 http:/ /pegasus. phast.umass.edu/ 

4 http://star-www.dur.ac.uk/ pdraper/gaia/gaia.html 

5 http:/ /www. eso.org/sci/observing/phase2/SMSpecial 
/RRMObservation.html 



2.10 KAIT 

The Katzman Automatic Imaging Tele scope (KAIT), also 
at the Lick Observatory |Li e t al. 2003), responded to the 
trigger and began taking unfiltered exposures starting 42 
seconds after the trigger time. This paper includes 206 un- 
filtered data points, which have been reduc ed in a standar d 
way and then calibrated to the i?-band (|Li et al.l 120031 ) . 
These data, once calibrated, are shown along with the 
PAIRITEL data in Figure H 

The first KAIT data point has a central time of 57 sec- 
onds, with a total exposure time of 30 seconds. Given the 
highly variable nature of early time GRB emission, the large 
error bars on the value, the long duration over which the 
magnitude was measured and (as later discussed) its outlier 
nature, we felt that this magnitude did not provide a useful 
measure the i?-band emission over this time. We therefore 
excluded it from the later analysis. 



http: / /www. eso.org/sci/software/esomidas/ 
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Figure 4. PAIRITEL and KAIT light curves for GRB 080310. 
The KAIT unfiltered data are calibrated to the ij-band. All data 
have been corrected for dust extinction. 



Table 1. ESO-VLT polarimctric observations of the late after- 
glow of GRB 080310 



t - t (s) 




Polarization (%) 


Band 


87171 


1447 


< 2.5 


V 


169501 


1447 


< 2.5 


V 


253724 


4607 


< 2.6 


V 



the physical be aming of the outflow jGhisellin i fc Lazzatil 
1999; Sari 1999). Therefore the polarization time-evolution 
is in principle a powerful diagnostic of the afterglow physics, 
and many attempts were carried out to comp are observa- 
tions to models (e.g.. lLazzati et aL 2003, 2004]) with partic- 
ular emphasis to the jet structure (|Rossi et al.ll2004T ). 

During the polarimetric observations of GRB 080310 
the afterglow showed a smooth decay (see Figure [3]) without 
any detectable temporal break or spectral change. Polariza- 
tion below ~ 2.5% cannot put specific constraints on the 
afterglow modeling or the jet structure. These results, are 
however compatible with what would be expected at late 
times, as this is when the forward shock should dominate 
emission and the magnetization signal of the fireball is lost 
in the interaction with the surrounding medium. 



2.11 Gemini 

Our last optical data were acquired with Gemini-North using 
the Gemini Multi-Object Spectrograph (GMOS) in imag- 
ing mode with the r-band filter. The observations began at 
10:22UT on 19th March 2008 and consisted of 5 x 150 second 
exposures. The data were reduced using the gemini-gmos 
routines within iraf. No significant flux was detected at the 
location of the afterglow, as reported in Table A5. 



2.12 Polarization 

Three linear polarimetric observation sets were carried out 
with the ESO-VLT during the late afterglow evolution as 
shown in Table[T] The observations allowed us to derive 
rather stringent upper limits althou gh still compatible w ith 
past late-time afterglow detection s |Covino et alj [20051 . see 
also, however, iBersier et al l 120031 ) and substa ntially lower 
than t he early time afterglow measurement bv lSteele et al.l 
|2009l ) for GRB 090102. In the case of GRB 090102, how- 
ever, the detectio n is taken at 160 seconds; a time when 
ISteele et all (|2009h argue that the flux should be dominated 
by reverse shock emission. 

The detection of linear polarization at the level of a 
few percent in the light from GRB optical afterglow is 
well within the predictio n of the external shock scenario 
l|Zhang fe M cszaros 2004j, and references therein) and in- 
deed it is still one of the most re levant observational findings 
supporting it (e.g.. lCovindl2010l 'l. On the other hand, a com- 
prehensive framework predicting the polarization degree and 
position angle evolution during the afterglow had to deal 
with the increasing complexity and variety of behaviours 
shown by the afterglow population. In general, the late- 
afterglow optical polarization is related to three main ingre- 
dients: the emission process able to generate highly polarized 
photons (i.e., synchrotron), the ultra-relativistic motion and 



2.13 Observations from literature 



Data from the Palomar 60 inch telescope (|Cenko et alfeoolf ) 
were obtained from the Palomar 60 inch- Swift Early Opti- 
cal Afterglow Catalog (ICenko et all 120091 '). in which the 29 
GRBs between the 1 st of April 2005 and the 31 st of March 
2008 with P60 observations beginning within t he first hour 
after t he initial Swift-BAT trigger are presented lCenko et al.l 
(2009) reduce data in th e iraf environment , using a cus- 
tom pipeline detailed in ICenko et al.l (|2006t ). Magnitudes 
were calculated using aperture photometry and calibration 
performed using the USNO-B1 catalog and the data were 
corrected for dust e xtinction using the extinction maps of 
jSchlegelet al.ll 19981 ). 

Further publi s hed d ata for GRB 080310 were obtained 
from iKann et all (|2010T ). in which SMARTS- ANDICAM 
data are detailed as part of an extensive survey of opti- 
cal data fo r GRB s in both the pie-Swift and Swift eras. 
IKann et al l <|2010T i reduce their data using standard pro- 
cedures in iraf and MIDAS. Both aperture and PSF pho- 
tometry were used in the derivation of magnitudes, when 
comparing to standard calibrator stars. 

The 0.6m Super-LOTIS (Livermore Optical Tran- 
sient Imaging System) telescope, lo cated at the Stew- 
ard O bservatory (Kitt Peak, Arizona; IPerez-Ramfrez et al.l 
|200J) began Ji-band observations of the error region of 
GRB 0803 10 at 08:38:43 UT 44 se conds after the start of 
the burst jMilne fc Williamsiliooi L The OT detected by 
IChornock et al.l ~ 20081 ) and confirmed by ICummings et all 
|200Sh was not apparent in the initial images, even when 
stacking the first three 10 second exposures. However, the 
subsequent 20 second exposures do show the optical tran- 
sient without stacking, which suggests that the GRB bright- 
ened in the f?-band during the first two minutes after de- 

7 http://www.nofs.navy.mil/data/fchpix 
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tection. A nearby USNO-B star was used to derive the R 
magnitude. 



3 MODELLING 

To model the emission of GRB 080310 we begin by fitting 
the BAT and XRT data, before then extending the model 
into the lower energy bands. 



3.1 Initial modelling of the high energy emission 

We ex pand on the previous work done by Willingal e et al.l 
(2010), where a sample of 12 GRBs were se l ected and fit- 
ted using the pulse model of lGenet fc Granotl (2009). In this 
model, the prompt emission component of GRB emission is 
split into a series of pulses, where each pulse is considered 
to be the result of a relativistically expanding thin spherical 
shell that emits isotropically. It was assu med that each o f 
the pulses was in the fast cooling regime (|Sari et al .11 19981 ) 
and that each of the X-ray and 7-ray spectra could be fit- 
ted w ith a temporally evolving Band function (Ban d et al.l 
Il993h . A Band spectral energy distribution is a smoothly 
broken power-law. Below, we show a modified Band func- 
tion, which includes temporal evolution. 



B{z) = B 



where 



E 



T-T e 



2 < 61 — hi, 
z > bi — 62, 



(1) 



(2) 



The times in Eqn [2] are all in the observed frame, where 
T e j is the time of shell ejection from the central engine and 
Tf is the time at which the last photon arrives from the 
shell along the line of sight. B is the normalisation to the 
Band function, with b\ and 62 being the low and high energy 
photon indices of the Band function, respectively. The spec- 
tral shape of each pulse evolves in time as det ermined by 
the temporal index d. In Willingal e et al.l (|20 10h d = -1 due 
to the assumption that the emission process is synchrotron 
under the standard internal shock model . It is this value 
of tem poral index that was also used by iGenet fc Granotl 
(2009) in the original derivation of the pulse profiles, mean- 
ing that only a value of d = -1 is strictly consistent with the 
original pulse model. 

The prompt pulse model describes the morphology of 
the prompt light curve and the rapid decay phase as ob- 
served in the high-energy bands. In addition to this, an af- 
terglow co mponent was also inclu ded in the modelling as 
outlined in IWillingale et~ai1 (|2007l ). The afterglow compo- 
nent has a functional form as outlined in Eqn [3] which com- 
prises an exponential phase that transitions into a power-law 
decay. 



/.(*) 



F a exp [Ota 



- T^J ex P 1 

ex p(^H 



t < Ta, 
t > T a . 



(3) 



In Eqn[3l f a (t) is the flux from the afterglow at time t, F a 
gives the flux at the transition time between the exponential 
and the power-law components, T a . T r is the rise time and 



Table 2. Properties of the prompt pulses identified in the BAT 
and XRT light curves. These include the peak time (T p j.), peak 
energy at this time (E p fS), the rise time (T r ) and arrival time of 



last photon (Tf) for each pulse. 


Pulse 


Tpk (s) 


E pk (kcV) 


T r (S) 


Tf 00 


61 


1 


-52.8 


200 


9.7 


44.5 


-1.49 


2 


-16.0 


200 


5.0 


6.0 


-1.20 


3 


-4.6 


200 


3.8 


11.7 


-0.40 


4 


1.8 


200 


2.7 


17.2 


-1.50 


5 


159.0 


12.3 


25.3 


74.0 


-0.30 


6 


191.6 


13.4 


12.4 


39.8 


-0.02 


7 


210.0 


21.3 


10.0 


46.7 


-0.16 


8 


235.0 


58.0 


8.0 


24.4 


-0.13 


9 


251.8 


42.0 


10.0 


52.4 


-0.20 


10 


282.0 


15.8 


10.0 


40.0 


-0.10 


11 


308.7 


16.1 


15.0 


32.5 


0.24 


12 


342.7 


7.8 


15.0 


40.1 


-0.18 


13 


366.0 


34.3 


10.0 


11.5 


-0.38 


14 


390.0 


1.2 


23.0 


56.1 


0.21 


15 


513.8 


3.8 


31.1 


193.8 


-1.26 


16 


582.1 


2.1 


43.9 


66.5 


-0.09 



finally a a is the index that governs the temporal decay of 
the power-law phase. 

Combining the prompt pulses and the afterglow compo- 
nent we adopted the same method of fittin g the data from 
the S wift BAT and XRT instruments as [Willingale et al.l 
(2010), by first identifying the individual pulses in the BAT 
light curve and allowing their parameters to be fitted by min- 
imizing the x 2 fit statistic for both the BAT and XRT light 
curves. The afterglow component was then fitted by allowing 
the routine to find the optimum values for the characteristic 
times and normalizing flux shown in Eqn [3] 

The simultaneous BAT and XRT fit is plotted in Fig- 
ure [S] and shows several important characteristics of both 
the burst and the model. Firstly, there is a lot of structure 
evident in the light curves, particularly during the first 1000 
seconds of the XRT light curves. In this fit, sixteen unique 
pulses have been identified and can be see n to fit the XRT 
data accurately. In the original fit from IWillingale et al.l 
(2010) there were only ten pulses, however this failed to 
fully model some of the structure in the softest BAT band 
(15-25 keV) between 200 and 400 seconds. The additional 
pulses now provide a better fit to this time, across all the 
BAT and XRT bands with a reduced \ 2 statistic of 1.68 for 
1212 degrees of freedom. The quality of this fit is formally 
not statistically acceptable, however the largest contribution 
to the x 2 value is from small scale intrinsic fluctuations in 
the data. The properties of the sixteen pulses are listed in 
Tabled 



3.2 Extrapolating to the optical and IR 

Previously, no optical or IR data have been included when 
modelling the light curve and a fit to only the Swift BAT and 
XRT data has been produced (Figure[S|. However, given the 
rapid response of optical and IR instruments to the trigger 
for GRB 080310 we extend the model to include these new 
sources of data. 

The simplest approach to fitting the optical and IR light 
curves for GRB 080310 was to use the fitting routine from 
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GRB080310 BAT 




GRB080310 XRT 




seconds since trigger 
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Figure 5. BAT and XRT data, showing the fitted I Willingale et ail fcoiol ) model. Left panels: BAT bands (Top to bottom: 100-350 keV, 
50-100 keV, 25-50 keV & 15-25 keV). Right panels: XRT bands (Top: XRT hard band, bottom: XRT soft band). Pulse parameters have 
been fixed, whilst the afterglow parameters have been fitted. Solid lines show the total fit (pulses and afterglow), whilst dashed lines 
show only the afterglow component to the fit. 



Willingale et al.l (|2010f ) and simply extrapolate the Band 
functions for both the pulses and also the afterglow to these 
lower energies. In this initial attempt all of the parameters 
previously discussed were held at the values obtained for the 
high energy fit, to see what modifications might be necessary 
to both components. 

Whilst it provides an acceptable fit to the XRT and 
BAT bands, the original pulse model vastly over predicts 
the optical and IR fluxes from the pulses (Figure [6]), which 
implies there must be a break in the pulse spectrum be- 



tween the X-ray and the optical and IR energies. Such a 
break is expected in a synchrotron spectrum, but can also 
been seen in thermal spectra as the Rayleigh- Jeans tail. Ad- 
ditionally, the afterglow prediction from the BAT and XRT 
fit is not consistent with the optical data at late times where 
the prompt component to the light curve is negligible. Figure 
[5]also shows that the temporal decay index of the power-law 
phase of the afterglow gives rise to a decay which is more 
rapid than the optical data suggest. 

An alternative method of reducing optical flux is to in- 
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Figure 6. Optical and IR light curves obtained by extending the high energy fit to the prompt pulses and afterglow component. The 
spectral shape of the pulses is assumed to be a Band function extrapolated back to the optical and IR data (K, H, J, I, R, V, B and U 
bands). Vertical dashed lines denote times of interest later discussed, at which the SEDs are considered at in more detail. The coloured 
dashed lines show the afterglow component to the extrapolation for each band. 



yoke d ust absorption, however, as reported in IPerlev et al.l 
(|2008j). extinction due to dust is low at A v = 0.10 ± 0.05 
at an average time of T (+l, 750) seconds for GRB 080310. 
As not only the late time emission seems unaffected by 
such absorption, but also by 2,000 seconds after trigger, we 
chose to favour a low energy spectral break in our mod- 
elling of the optical emission. The following sections of this 
paper detail the imple mentation of such alterations to the 
IWillingale et all ||201Ch model. 

3.3 Modifications to the model 

3.3.1 An additional spectral break 

To reduce the flux from each pulse in the optical and IR 
light curves of GRB 080310 we introduced an additional 
break to the spectrum for each prompt pulse. An example 
of such a spectrum is shown in Figure [7] where a regular 
Band function has a low energy break below E pk . To fully 
describe this break, three parameters are needed; a value of 
the break energy (E a ), the spectral index of the power-law 
slope in the low energy regime (p a ) and a temporal index 
which describes how the break energy evolves in time (d a )- 
The value of E a is defined at T p k, when the emission from 



the pulse is at a maximum. The entire pulse spectrum al- 
ready evolves in time, according to the index d as shown 
in Eqn [2l and so the expectation was that the time evolu- 
tion of E a would be related. Previously a value of d = -1 was 
used, which when integrating of the Equal Arrival Time Sur- 
face (EATS) means a pure Band function can be recovered 
in the source frame given a Band function in the observed 
frame. Other values of d prevent assumptions to be made 
concerning the source spectrum. We allowed the index d a to 
be fitted independently, to test whether the evolution of the 
break was the same as the rest of the spectrum. Again, we 
note that only d a = -1 is fully consistent with the original 
I Genet fc GranotJ (2009) pulse model. The functional form of 
the new spectral model for each pulse is shown in Eqn 2] 

( AE Pa ~ 1 exp E~ a , E < E a , 

N(E) = I AE^- 1 exp (^f) , E a < E < E pk , 

{ AE" 2 " 1 exp (a - 0) E^ a - p) , E > E pk . 

The model spectrum for photon emission assumes several 
things. Firstly, there is a single underlying population of 
relativistic electrons, whose energies can be described by a 
broken power law such as that described in Sh en fc Z hane 
(2009). Such a spectrum of electron energies leads to a simi- 
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Figure 7. Example modified spectrum including additional spec- 
tral break. 

lar photon spectrum: a singly broken power-law, which cor- 
responds to the observed Band function. However, it would 
be unphysical for this spectral shape to extend indefinitely 
to low energies, particularly as the electron energy spectrum 
has an associated minimum energy E m . This minimum elec- 
tron energy has a related emission energy at a frequency v m . 
With no self-absorption, the spectral index to the emission 
spectrum changes to | at photon frequencies below v m , but 
if self-absorption is present then at low energies we would 
expect to see a steeper spectral index of 2 if the absorption 
frequency is less than v m or alternatively a spectral index 
of | in the case of v m < v a . In the former case, if there is 
an intervening spectral range between v a and v m , a spectral 
index of | is expected between these two frequencies. 

To keep the model simple, due to the limited nature of 
the spectral information available, we only included a single 
additional break in the spectrum. By allowing the resultant 
spectral index to be constrained by the observations, rather 
than assuming a value, we hope to better understand the 
physics required to explain the spectrum of prompt GRB 
pulses. 

3.3.2 The afterglow 

As seen in Figure [Bl prompt pulses without a low energy 
break significantly over predict the observed flux in the op- 
tical and IR light curves prior to 10,000 seconds. By intro- 
ducing a low energy break to the pulse spectra, we hoped 
to completely remove the prompt contribution to the total 
emission, and model the entire light curve with afterglow 
emission. One reason for considering this is the smooth na- 
ture of the optical light curves. Whilst there is some varia- 
tion during the plateau seen between 200 and 3,000 seconds, 
this is at a low level, and happens smoothly over a large pe- 
riod of time. If the prompt pulses observed at higher energies 
were also the dominant source of emission at these energies, 
then we might expect to see similar structure in the optical 
and IR light curves to that shown in Figure [5] 

The first issue to address with the afterglow fit was 
to account for the small difference in temporal decay in- 
dex (a a ) between the higher and lower energy bands. To 
do so, we included an additional parameter that described 
the difference in this index, and allowed it to be fitted. By 
including this, we could account for the slightly shallower 
decay of the power-law phase in the optical and IR chan- 



nels. The change in a a required was found to be small, at a 
value of approximately 0.2, but the associated errors in each 
instance showed it to be inconsistent with zero. This may be 
explained physically by a slight curvature in the spectrum 
of the afterglow between the optical and X-ray bands. This, 
however, is not a large enough difference to suggest a spec- 
tral break, such as the one introduced to the prompt pulses. 
The GRB afterglow flux is also thought to be synchrotron 
emission. However, at these late times the radiation is usu- 
ally assumed to come from an optically thin plasma, and so 
a self-absorption break would be expected at energies lower 
than the optical bands observed for GRB 080310. We also 
ruled out contamination in the optical and IR wavebands as 
being the source of this difference as any emission from the 
host should be constant with time. This should add a con- 
stant offset to the data, and should the afterglow emission 
reach a comparable order of magnitude, a plateau at the 
end of the optical decay would be expected. This plateau is 
not observed, and the difference between the low and high 
energy temporal decay indices of the afterglow is also deter- 
mined by data prior to when the host would be seen to make 
a significant impact on the optical and IR light curves. 

With no prompt component capable of rising quickly, 
we also had to modify the manner in which the afterglow 
rises, to account for the rapid increase in flux seen in the 
V-band data at approximately 100 seconds in Figure [3] To 
do this, we introduced a third part to the functional form 
that is shown Eqn [3] This extra regime is shown in the top 
line of Eqn [5] and describes a power-law rise in flux, with 
an index of 2. We also allowed the afterglow to be launched 
at a time that was independent to the trigger time of the 
GRB, by introducing a launch time (T;), which offsets the 
afterglow in time to better fit the timing of the rise observed 
in the V-band. These modifications are shown in Eqn [5] and 
are used for the afterglow-dominated modelling only. In the 
case of t he prompt dominated fit , we return to the afterglow 
model of IWillingale et alj (|2007T ) as explained in §3.5. 

{F a exp (a a - ^f^" ) (w^k) * > t < T r , 
F a exp (a a - (i ~ a -T"° ) , T r <t<T a , 

Fa ( Ta ^ , t > T a . 

3.4 Afterglow dominated fit 

Given the smooth nature of the optical and IR light curves, 
our initial use of the low-energy spectral break in the prompt 
pulse spectrum was to remove the prompt component to the 
optical light curves entirely, and try and explain the early 
time behaviour with only afterglow emission. To do this, the 
prompt pulses were initially switched off from the optical 
and near infrared fit. With these removed, the afterglow pa- 
rameters T(, T r , T a , F a , OL a and the change between optical 
and high energy values of a a were fitted. Having obtained 
values for these parameters using the fitting routine, we re- 
introduced the prompt pulses, allowing p a to be fitted, to 
find the minimum break in the pulse spectra required to 
remove the prompt optical and IR flux. Where the fit statis- 
tic tended to unphysical values of p a we checked the effects 
of forcing p a to be | . When this was undertaken, we found 
that whilst there was an increase in the reduced \ 2 statistic, 
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this was incredibly small, approximately 1 x 10 3 ,for 1642 
degrees of freedom. 

3.4-1 Light Curves 

Figure [8] shows the light curves for the best fit found us- 
ing the method described above, for which a reduced \ 2 
value of 2.44 with 1,642 degrees of freedom was found. The 
late time XRT afterglow is n't as well fitted as in the orig- 
inal IwniingalejjElil ||201(t l method (Figure [5]), where only 
the high energy bands were included. This is particularly 
the case for the last data points, and those at the tran- 
sition between the prompt and the afterglow phase in the 
X-ray bands. The cause of this is due to the greater number 
of optical points being the largest constraint on the after- 
glow parameters, particularly the characteristic times, T a 
and T r . Because these times have changed from the origi- 
nal fit, the temporal decay index in the power- law phases 
has reduced from 1.58±^ to 1.30±g . The afterglow pa- 
rameters all have a reduced error in this newer fit as the 
optical and IR data points allow for more accurate fitting 
with a greater quantity of data. The characteristic times 
also have improved accuracy, with T r = 223^20 seconds, T a 
= 6, 368lgQ 5 seconds and Ti — 1201^ seconds. The errors 
quoted in this instance were obtained when floating all of 
the parameters mentioned simultaneously. Once the after- 
glow parameters were in place, the final bulk value of p a 
was found to tend to 5.6, with a lcr lower limit of 5.5 and 
an unconstrained upper limit. The lower limit is tightly con- 
straining, and likely due to the fitting routine trying to re- 
move contributions from the brightest early pulses, which 
therefore affects the value assigned to the weaker pulses as 
well. The lack of an upper limit suggests that the fitting rou- 
tine is attempting to remove the prompt component entirely 
from the optical emission, as driven by the morphology of 
the light curve particularly the dense sample of KAIT data 
points. The value of 5.6 is very significantly steeper than that 
expected at either 2 or | . The suggested lower limit from the 
X 2 distribution is surprisingly tight, given the shape of the 
curve obtained for the statistic when only varying p a as the 
distribution is remarkably flat for a broad range of values. 
Looking at the fit statistic by eye, it would appear that the 
expected values from synchrotron or thermal spectra are not 
entirely ruled out, but lie within the lower limits of where 
the fit is a faithful representation of the data. To test this, 
we forced p a to a value of | for all the prompt pulses, and 
noticed that the value for the \ 2 statistic increased by a 
nearly negligible amount. By using p a = § self-absorption 
becomes a viable explanation for the physical mechanism 
causing the observed level of optical and near infrared flux. 
It is therefore this value that was adopted for p a for the 
prompt pulses when trying to suppress them at early times 
in the lower energy regimes. 

The broad scale morphology of the optical and NIR 
light curves are described by the fast rising afterglow, with 
the rise observed in the U-band being fitted and the level of 
flux of the plateau seen between 100 and 1,000 seconds being 
consistently modelled in five of the six bands in which it can 
be seen. The exception to this is that the [/-band model 
light curve morphology seems to over-predict the emission 
actually observed in the first 1,000 seconds, which may be 
a result of the correction required to removed the effects 



of absorption from the Lyman forest as described in the 
section detailing the UVOT observations. The I-, R- and 
1^-band data during this plateau phase do appear to have a 
slight dip in flux, which is not picked up by the afterglow 
fit. Additionally, the first V-band points indicate a decline 
in flux prior to the launching of the afterglow. To fit this, 
one of the early time pulses would have to be switched on at 
low energies, implying that at the earliest times the prompt 
component is still important in the optical and IR regimes. 

The V-band data between 100 and 1,000 seconds are 
from three sources, VLT V-band, UVOT v-band and the 
UVOT white filter, which is normalised to the «-band. This 
normalisation is done at a time when there is near simul- 
taneous UVOT v-band and white filter coverage. The white 
data are then adjusted at this time, to make the observed 
flux the same as that in the v-band. The same correction fac- 
tor is then applied to all the white filter points, which relies 
on an intrinsic assumption that the spectrum is not varying 
between the u-band and white filter for all observations. If 
the total GRB low energy emission is primarily from the af- 
terglow, this assumption is accurate, as the spectral break 
in the Band function will not migrate to the optical part of 
the spectrum. We noted, upon inspection of the data, that 
the UVOT points between 300 and 2,000 seconds appear to 
not be entirely consistent with the near simultaneous VLT 
data. The VLT flux values have much lower associated er- 
ror, which led to us reprocess the UVOT data, to verify the 
values. Doing this we obtained the same UVOT fluxes. To 
account for the discrepancy, we introduced a systematic er- 
ror in the UVOT points at the level of 10% of the measured 
flux for each datapoint. Analysis of only the VLT data, still 
shows the deficit of flux, which we discuss throughout this 
work, at a level that is significant given the small errors 
associated with the data. As this feature is significant in 
instruments other than UVOT, we believe that the discrep- 
ancy from the UVOT instrument does not impinge on the 
main conclusions of this work. 

Figure [3] would suggest that the decline in flux is not 
marked by a change of instrument. 

The temporal index d a governing the time evolution of 
E a is tightly constrained in this model, to a value of d a = 
— 1.011q'q3- Not only does this mat ch the temporal evolutio n 
of the characteristic energy E c in iGenet fc GranotJ {2009), 
but the 1-cr errors quoted are restricting. Having analysed 
the one dimensional \ 2 surface, it is clear that a more rapid 
evolution of the temporal index allows the prompt pulses to 
peak at earlier times, when the normalising flux at the peak 
of each pulse is still significant, causing an over-prediction of 
early time flux. At lower values the \ 2 distribution is fairly 
flat, as the prompt flux normalisation is so low, it cannot be 
seen to contribute to the total light curve. 

3.4-2 Spectral Energy Distributions 

Another way of comparing the model to the observed data 
is to consider the spectral energy distribution (SED) at fidu- 
cial times. The times selected are indicated in Figures [5] [5] 
and [10] as dashed vertical lines. These correspond to times 
of 350 seconds (where there are early-time optical data and 
the activity in both the X-ray and 7-ray bands), 1,800 sec- 
onds (a time where the high energy only model shown in 
Figure [6] suggests the afterglow is becoming significant), 2.5 
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Figure 8. Simultaneous fit to optical, IR and XRT data, in which the early time optical and IR are fitted by the afterglow emission. 
Top panel: Optical and IR bands fitted with low energy break (K, H, J, /, R, V, B and U bands). Bottom left: XRT soft band fit. 
Bottom right: XRT hard band fit. The afterglow parameters are fitted, as are values for d a , whilst E a = 0.3 keV and p a has been set to 
|. Coloured dashed lines show the afterglow component in each band. The reduced \ 2 for this fit is 2.44 with 1,642 degrees of freedom. 



x 10 4 and 1 x 10 5 seconds (both of which are later times 
that should be dominated by the afterglow). Figure [5] shows 
the model SED at these times, and also those at both 200 
and 1,000 seconds. The times shown were chosen to both 
allow for comparison with the maximum number of bands 
(both optical and X-ray) at each time and to show times 
where the SED can discriminate between prompt and after- 
glow emission at optical and IR energies. The SEDs shown in 
Figure [9] are those for the fit where the afterglow dominates 
the entire low energy emission. In each case, error bars are 
plotted, but in some panels are too small to be seen. 

Figure [9] shows that by combining the prompt pulses 
and the afterglow component, the model SEDs have compli- 
cated structure at high energies. The first panel shows the 
SED at 200 seconds, a time by which the first six pulses 
have peaked. As an example of time evolution, the spectral 



contribution from pulse four is shown with the dashed line 
in this and subsequent panels. The temporal evolution of 
pulses 4, 12 and 15 can be traced through all six panels. At 
200 seconds pulse 4 is the one of the three pulses causing de- 
viation from the afterglow Band function (the others being 
pulses 1 and 6). Following it through all six panels of Figure 
[9]shows the general evolution of all the pulses. The peak flux 
reduces with time, and the energy at which this peak is seen 
becomes lower with increasing time also. By the late pan- 
els, the component of pulse 4 can still be seen on the axes, 
however, being several orders of magnitude lower than the 
afterglow Band function, it does not produce an observable 
feature in the total SED. Because the spectral index for all 
the pulses is at a value of p a — 2.5, before the peaks migrate 
to the optical regime, their contribution to the total flux 
at these energies is negligible. Once the peak has had time 
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Figure 9. SEDs for fiducial times in the light curve of GRB 080310 using the afterglow model to fit the early time, low energy emission. 
Prompt pulses are included with a bulk value of p a = 2.5. The solid line is the total SED and the dotted line is the contribution from 
the afterglow component. Pulses 4, 12 and 15 are also shown by the dashed and dotted lines as outlined in the key shown in the bottom 
right panel. These are included as examples to show the evolution of some of the longer lasting pulses in the fit. The colours of the data 
points follow the same scheme as that used in the light curves, with the BAT and XRT data being represented by lines which extend 
along their entire coverage. 
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to evolve to such energies, the normalizing flux has been 
diminished, so the pulse contributions remain negligible. 

Figure [9] also shows that the dominant pulses in the 
SED of GRB 080310 do not have to be those most recently 
launched. Pulse 12 is emitted after pulse 4, as shown in the 
first two panels of the figure. Despite this, pulse 12 has dis- 
appeared entirely by 2.5 x 10 4 seconds, whilst the earlier 
pulse can still be seen on the axes, even though the mod- 
elled emission is at a level significantly below the total SED. 
The bottom two panels show late times in the light curve, 
where only the afterglow contributes flux to the total SED. 
Looking only at the afterglow in all six panels, it can be seen 
that it varies very little at early times. This is due to the 
combination of a quick rise time after an early launch with a 
late transition time between the exponential and power-law 
decay phase for the afterglow component. This implies that 
there is a period where the afterglow effectively plateaus for 
a few thousand seconds. This behaviour can be identified in 
the light curves shown in Figure [H] 

The SEDs in Figure [9] suggest that the optical emis- 
sion is fitted by the extrapolation of the afterglow Band 
function, particularly at late times. The second panel is per- 
haps the most puzzling as the optical and near infrared data 
are aligned approximately with the afterglow SED, however 
there is a clear dip in the light curves, which suggests that 
a component which simply rises then falls smoothly can- 
not explain the morphology that is seen. We have not tried 
to fit a low energy break to the afterglow Band function, 
as discussed previously, because such a break is expected 
at energies well below the //-band, due to the optically thin 
nature of the circumburst medium. Additionally, the correct 
level of flux has been reproduced without the introduction 
of a further spectral break. 

3.5 Early time prompt dominated fit 

By fitting the optical and IR emission with only a significant 
afterglow component, there were discrepancies between the 
model and the data. Namely the potential decay before the 
sharp rise observed in the V-band and the deficit in flux ob- 
served in the V- and _R-bands between 300 and 2,000 sec- 
onds. In addition to this, the ?7-band data are over predicted, 
which is possibly due to the absorption correction from the 
Lyman forest as previously discussed, or alternatively could 
be better described with an alternative model. In an attempt 
to explain these features, we considered a model where these 
early-times were dominated by the prompt emission in the 
optical and near infrared regimes. 

3.5.1 Light Curves 

We first only considered the prompt pulses, by turning off 
the afterglow component and fitting only prior to 2,000 sec- 
onds. Once the early time emission was fitted, we then in- 
cluded the afterglow to account for the late time emission. 
The potential for degeneracy between p a and E a was con- 
sidered, as the energy difference between the X-ray data and 
the optical and IR bands is over two orders of magnitude, 
while the U- and Jf-bands are only separated by less than 
a single order of magnitude. Data were not taken simulta- 
neously in these two bands which are the most spectrally 



separated. It is therefore difficult to determine both E a and 
p a independently. For this reason the break energy was set to 
a value of 0.3 keV (at T p k) which is the soft end of the XRT 
spectral range. This implies that the fit to the high-energy 
regime remains unaltered, whilst giving the maximum en- 
ergy range in the SED between E a and the optical and IR 
bands. Therefore the value for p a obtained is a lower limit, 
as moving E a to lower values would steepen the power-law 
index in the low energy part of the spectrum. 

Initial attempts to model the low energy flux with 
prompt pulses held the parameters p a , d a and E a to the 
same value for all pulses. This led to a bulk value of p a = 
0.65^qq2> which is significantly different from the spectral 
index of the Band function below E p k, but not consistent 
with the expected values of p a if self-absorption was having 
an effect on the emission that is observed. The ler limits on 
this value do not include |, which is what would be expected 
if there is an absence of self-absorption, and the optical and 
IR bands are below the emission frequency of electrons with 
the minimum energy in the electron population. The reduced 
X 2 value is 3.61 in this instance, with 1,642 degrees of free- 
dom. There were still some issues with this fit, such as the 
fit not describing the slight decay followed by sharp rise in 
flux seen in the V-band. So we decided to allow all pulses to 
have independent values of p a , to see if this could account 
for the features seen in the data. The method adopted to do 
this remained the same as that for the bulk value of p a , and 
the light curves obtained are shown in Figure [lOl 

Given the ch oice of two afterg l ow m odels; the original 
functional form of IWillingale et all (|2007T ) and the modified 
version shown in Eqn[S] we found that in this instance, the 
better fit was obtained with the former, as the smoother 
peak to the afterglow component allowed for a closer fit to 
some of the J-, K- and //-band data at this time. This also 
removed a free parameter from the fitting routine. We still 
retained the difference between optical and high energy a a 
to describe the late time decay. With these modifications, 
the early time prompt dominated model had a reduced \ 2 
of 2.43 for 1,642 degrees of freedom. This is significantly 
better than the prompt model with a bulk value of p a for 
all pulses, with a change in total \ 2 °f over 1,850. 

There are little data to constrain p a for the first four 
pulses, which are predominantly observed with BAT, and 
so these values have large associated uncertainties. Despite 
this, it was clear that a large value of p a was required for 
pulses 1, 2 and 4, to allow for the decay and subsequent rise 
seen in the V-band data. The best fit to the initial decay 
was found by modelling it with pulse two, and suppressing 
pulses one and three so that they did not contribute to the 
optical and IR emission observed from the beginning of the 
optical coverage. When actually fitted, the values of p a for 
pulses 1, 2 and 4 all were larger than |, and so in the same 
manner as all the prompt pulses in the afterglow-dominated 
model, these were set to p a = | as Figure [TO] shows that 
the sharp rise is now well fitted by the rise of pulse four 
in a similar way to the simultaneous rise seen in the XRT 
data. The data obtained from KAIT have an earlier point 
at 57 seconds, which does not show this decay, it is difficult 
to fit this individual point, as whilst it is at an ideal time 
to investigate the prompt optical behaviour of the GRB, it 
is a 30 second observation, and so it is an average over a 
time during which large variations might be expected. This 
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Figure 10. Simultaneous fit to optical, IR and XRT data, in which the early time optical and IR are fitted by the prompt emission. Top 
panel: Optical and IR bands fitted with low energy break (K, H, J, I, R, V, B and U bands). Bottom left: XRT soft band fit. Bottom 
right: XRT hard band fit. E a = 0.3 keV, p a and d a are allowed to float (producing a value of d a = -1.01) and the afterglow parameters 
(T r , T a , oi a and F a ) are also fitted. The reduced \ 2 for this fit is 2.43 with 1,642 de grees of freedom. 



is why we have not tried to fully explain this point with 
our model and have removed it from the fitting undertaken. 
Indee d, non-detections from both Super-LOTIS and RAP- 
TOR (jWozniak et alJl200St ) at similar times to this KAIT 
datapoint offer further evidence towards the faintness of the 
emission in the optical regime at this time. 

Analysis of the residuals for the optical, IR, XRT and 
BAT data show there are three clear contributions to the 
total x statistic. The first contribution to the fit statistic 
is from the XRT afterglow component of the model. This 
deviates from the observed data at two times. Firstly, the 
transition from the rapid decay phase to the afterglow shows 
an over prediction by the model, as do the data after 10 5 
seconds. The latter of the two features could be explained 
by a late temporal break in the X-ray afterglow. Such a 
break should be achromatic, and therefore also observed in 



the late-time optical data. The final _R-band observation, at 
approximately nine days after the trigger, is fainter than 
the extrapolation from the final power-law decay in the 
afterglow-dominated model, which would support the pos- 
sibility of a jet break. With the prompt-dominated optical 
and IR model, there is a suggestion of either a plateau or 
rebrightening at approximately one day after the trigger, 
followed by a steeper decay. In both models, the data are 
too sparse to justify the inclusion of further parameters in 
the late-time modelling of the optical data. 

The next significant contribution to the value of \ 2 
comes from the early-time REM data points between 150 
and 600 seconds in the if-band. These are under predicted 
by the model as shown in Figure fTOl Due to the sparse nature 
of the NIR data at this time, the fit is driven by the more 
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densely sampled V- and -R-band data. This discrepancy is 
better highlighted in the SEDs shown in Figure [TT1 

The other notable source of contribution to the fit 
statistic is the KAIT data. The data were initially reported 
in unfiltered magnitudes, with values ranging between 17 m 
and 20 m , and errors typically less than a tenth of a mag- 
nitude. With such small errors for over two hundred data 
points the contribution to the total validity of the fit from 
the intrinsic scatter of the dataset is significant, despite the 
fact that the data and model can be seen to agree by eye. 



3.5.2 Spectral energy distributions 

As can be seen in the SEDs shown in Fig. 1111 by introduc- 
ing the prompt component to the spectrum of GRB 080310, 
each SED has structure beyond the simple Band function of 
the afterglow at earlier times in the low energy bands. For 
the prompt-dominated model of the early-time low energy 
emission, the spectral indices typically lie in the range -0.65 
< p a < 1- This does not include the earliest pulses, which 
have a necessarily steep spectral index to prevent their emis- 
sion dominating (and poorly fitting) the optical and IR light 
curves. These have been fixed to a value of p a = § , as well 
as pulses 8 and 16, which also require steep indices. 

Unlike in the afterglow-dominated model, the SED at 
350 seconds shows the optical and IR (particularly the H- 
band) model to disagree with the data; this is the most valu- 
able of the SEDs in Figure 1111 as it is one of the only two 
prior to the afterglow making the dominant contribution 
to the total low energy emission, and of these two has the 
best low energy coverage. The R-, I- and V-band data at 
this time are at approximately the correct level of flux for 
each band, although the /-band and _R-band deviate from 
the expected shape of the SED. However the //-band data 
(which offers the largest range of spectral information when 
considered with the other three bands) is noticeably under 
predicted, as corroborated by the light curves in Figure [TDl 
When inspected by eye, the four bands at this time could 
be thought to lie on a single power-law with a spectral in- 
dex similar to that of the afterglow. This would lend more 
credence to the afterglow -dominated early-time optical and 
near infrared model. 

The subsequent SEDs are all dominated by the Band 
function of the afterglow at optical and IR energie s, and, 
despite the use of the orig inalfWillingale et alj |200it ) after- 
glow model and changes in T a , F a , and T r , the late time fit 
is similar to that obtained for the afterglow -dominated fit. 
The optical and near infrared data points are all over pre- 
dicted at 2.5 x 10 5 seconds, which something that can be 
seen by looking at the light curves, particularly in the case of 
the UVOT u-band data. The it-band dataset displays some 
unusual characteristics, being significantly fainter than ex- 
pected prior to 1,000 seconds and then appearing to plateau 
between 2 x 10 4 and 1 x 10 5 seconds. We considered the 
removal of this dataset from the analysis, but as both the 
values of parameters obtained for the best fit seem largely 
insensitive to the it-band data we retained them. This is a 
result of the relatively sparse coverage in the it-band, prin- 
cipally in comparison to that of the /?-band, the weighting 
of the u-band data in determining the fit statistic is small. 



4 DISCUSSION 

For direct comparison between the two models discussed in 
the previous section, their parameters have been included in 
Table [3] Given the degeneracy between E a and p a it isn't 
possible to simultaneously fit exact values to both without 
a more extensive low energy coverage over a larger range of 
wavelengths. Because of this, E a was fixed at 0.3 keV, in an 
attempt to find a value of p a . 

There is no clear difference in the fit statistic for either 
of the models presented in this work. The difference in to- 
tal x 2 is l ess than 30 over such a large range of degrees of 
freedom. Whilst numerically, this difference is in favour of 
the prompt-dominated early-time optical and near infrared 
model, we do not believe that the magnitude of the differ- 
ence is enough to warrant favouring one model above the 
other. 

For the afterglow-dominated early-time model the fit- 
ted bulk low-energy spectral index for the prompt pulses is 
unphysically steep (p a = 5.6). However, the change in fit 
statistic when forcing this value to one consistent with self- 
absorption was near negligible, which is why we therefore 
adopted the less extreme value. In contrast, most values of 
p a are more reasonable (-0.65 < p a < 1) when the early-time 
low-energy emission is dominated by the prompt pulses iden- 
tified at higher energies. This is not, however, true for all the 
pulses. It was necessary to suppress pulses 1, 2, 4, 8 and 16 
with positive, steep low energy spectral indices so that the 
optical and IR emission was not over predicted. To under- 
stand this behaviour we looked at the one dimensional x 2 
distribution for p a of the earliest of these three pulses. We 
found that the distribution reduced asymptotically at higher 
values of p a , meaning that there was no clear minimum. 
However, the actual reduction in the fit statistic achieved 
by reducing the low energy spectral index rapidly dropped, 
so this decrease in \ 2 value quickly stopped making a sig- 
nificant difference to the quality of the total fit compared 
to other contributions to the \ 2 statistic. After establishing 
this, we again adopted values of p a = | in these instances, 
with a minimal increase of \ 2 being introduced as a result. 

Given the variations in p a for the prompt dominated 
model, we tried to relate these values to other parameters 
of the pulses, such as pulse duration (T p k) and the low en- 
ergy photon index of the Band functions for each pulse (6i). 
We found no correlations between p a and the other pulse 
parameters. 

In the afterglow dominated fit, the prompt component is 
several orders of magnitude fainter than the observed emis- 
sion. Additionally, with positive values of p a , time is required 
for the peak of the pulse spectrum to migrate to the opti- 
cal and IR bands, before which the spectrum is rising at 
these energies. The pulse emission is largest around E a , the 
evolution of which is governed by the temporal index d a , 
which has been shown to be negative. This means that the 
pulse energy decreases with increasing time. In addition to 
this, the normalization of each pulse scales with t , which 
means that when the peak of emission reaches the optical 
part of the spectrum the total flux is reduced. To illustrate 
the prompt component in this scenario, we subtracted the 
afterglow component from the afterglow-dominated model, 
and produced Figure [P2l which shows the //-band prompt- 
only light curve that underlies the afterglow emission. 
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Figure 11. SEDs for fiducial times in the light curve of GRB 080310 using the prompt-dominated early time model. The solid line is 
the total SED and the dotted line is the contribution from the afterglow component. Pulses 4, 12 and 15 are also shown by the dashed 
and dotted lines as outlined in the key shown in the bottom right panel. These are included as examples to show the evolution of some 
of the longer lasting pulses in the fit. The colours of the data points follow the same scheme as that used in the light curves, with the 
BAT and XRT data being represented by lines which extend along their entire coverage. 
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Table 3. Pulse p a values and afterglow parameters for both the prompt and afterglow dominated ea rly time optical emission models. 
Note that the value for T; in the prompt dominated model is untabulated, as the Willinealc ct al. ( 2007) model was used for the afterglow 
in this instance. All instances where p a = | have the values fixed rather than fitted, and so no uncertainties were calculated for them. 
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Figure 12. Modelled i/-band light curve showing the prompt 
component only for the afterglow dominated fit. 



The afterglow-dominated, early-time optical model re- 
quired three alterations to the afterglow component; a 
change of temporal power-law index between high and low 
energies, a power-law rise at early times and a variable 
launch time T;. However, even with these modifications, 
there were several features of the data that were not entirely 
modelled in this scenario, including the apparent decay in 
the V-band prior to the rapid rise and the slight deficit in 
flux between 300 and 2,000 seconds. The tendency towards 
unphysical values of p a also suggested that fitting the early- 



time emission with a fast rising afterglow may not provide 
the best model of the observed flux. The SEDs of Figure 
[5] confirm that the afterglow-dominated fit provides a rea- 
sonable fit to the data and that at 350 seconds, which is 
a time where at high energies the prompt components are 
dominant, the afterglow Band function fits both the level 
and slope of the SED. 

When fitting the early-time flux with the prompt pulse 
emission model we first adopted a single value of p a which 
was held constant for all pulses. The low energy spectral in- 
dex obtained was p a =0.65, however most of the early time 
emission was hidden by the tail of one of the earliest prompt 
pulses. Aside from a poor fit to the ff-band data, the fast rise 
seen in the l/-band was not well represented in this instance, 
so the pulses were allowed to have independent values of p a - 
This produced a better fit to the data, significantly improv- 
ing the x 2 statistic. The fast rise in the F-band, and the 
small scale variations in the light curves seem to be rep- 
resented well by this model. There were, however, discrep- 
ancies with this alternative. Firstly, the //-band data were 
under predicted by the model, which can be seen in both 
the light curves and SEDs shown for this model. The SEDs 
also show that, at 350 seconds, the optical SED appears to 
be more accurately represented by a single power-law which 
agrees more with the nature of the afterglow. Three of the 
first four pulses (and also pulses 8 and 16) tend towards an 
unphysically steep spectral index when fitted. As with the 
afterglow-dominated model, constraining the values of p a to 
I in these instances does not significantly reduce the quality 
of the fit obtained. 

The other values of p a for the prompt-dominated fit 
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could be consistent with the energy bands lying below the 
energy of emission from the least energetic electron in the 
relativistic population responsible for the observed photons. 
We found no relation between the low energy spectral in- 
dex for each pulse and any other pulse parameter, giving 
no insight into the cause of the variation of p a - An inter- 
esting result is the marked difference between pulses 1, 2, 
4, 8 and 16 when compared with all the other pulses. p a 
for these early pulses indicates their emission is from an 
optically thick environment in the optical regime. Another 
implication is that these pulses must be launched from a 
different environment than the other eleven which have far 
shallower spectral indices at low energies. 

One additional benefit to the prompt-dominated model 
was the ability to retu rn to the simpler afterglow model of 
IWillingale etafl (|2007r i, as the rapid rise seen in the V-band 
was attributed to the rise of pulse four. 

An alteration attempted with the prompt-dominated 
early-time model was to introduce an offset to the onset of 
the afterglow. It was hoped that by allowing the afterglow to 
rise earlier, the SED at 350 seconds could reconciled with the 
model. This was found not to be possible, as this returned 
the model to something resembling the afterglow-dominated 
fit, and therefore we were unable to model the small scale 
variations and observed optical deficit in the plateau be- 
tween 300 and 2,000 seconds. 

Both of the alternative models presented in Figures [5] 
to [11] are poor at fitting the transition between the rapid 
decay phase and the afterglow plateau observed in the X- 
rays. We attempted to better model this by allowing the 
afterglow to rise at a later time than presented in either 
of the two fits highlighted in Table [21 By doing so, it was 
possible to improve the quality of the high energy afterglow 
fit, but to the detriment of that obtained for the optical 
data. With an afterglow that rises later but more quickly, 
the optical data between 3,000 seconds and 8,000 seconds 
are under predicted by the model suggested. With a larger 
number of data points, and correspondingly better statis- 
tics, the optical points were those that we therefore favoured. 
The prompt-dominated, early-time optical and IR model de- 
scribed in Tableland illustrated in Figures [TO] and [TT] sum- 
marize this model. Unfortunately, the XRT coverage con- 
tains a gap during the afterglow plateau phase, and therefore 
the exact morphology of this component, which would help 
significantly constrain the rise of the afterglow, is unknown 
at these times. 

An alternative suggestion, given the acceptable fit at 
approximately 2,000 seconds in all the available optical IR 
bands, to explain the factor of two or three difference be- 
tween the X-ray data and model at this time would be to 
include a spectral break in the afterglow spectrum. Whilst 
this may improve the fit at these times, the later afterglow 
coverage between 10 4 and 10 s seconds show both the higher 
and lower energy bands to be modelled at the correct level. 

Additionally, a large contribution to the \ 2 distribu- 
tion was from the optical and IR data, particularly the 
KAIT dataset, with low magnitude errors. Whilst the ob- 
servational errors in these datapoints may be as reported, 
each of the datasets had to be calibrated so they all were 
in standard bands. In doing so the systematic errors associ- 
ated with these data points increased. We considered this, 
and added a systematic of 0.03 m to every optical or IR point. 



Table 4. t peak values and initial bulk Lorentz factors for both 
the prompt and afterglow- dominated early-time optical emission 
models in both an ISM and wind like environment. 
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The quoted values for the fit statistics include this system- 
atic source of error. 

Knowing the characteristic times of the afterglow for 
each of the two models allows us to calculate the initial bulk 
Lorentz factor in both an interstellar medium (ISM) like or 
wind do minated circumburst environ ment. To do so we u sed 
Eqns[6l (|Molinari et al.ll2007T ) and[7l (|Sari fc Piradll999h . 



To, 



/ 3E lso (l + z) 3 
\32TTnm p c 5 V t 3 peak 

( E lso (1 + z) 
y8TvAm p c 3 r]t p ,, ak 



(6) 



(7) 



In the equations above Ei SO is the isotropic equivalent en- 
ergy of the GRB, r\ is the radiative efficiency of the fire- 
ball, z is the measured redshift, n is the number density of 
the circumburst medium and A is a normalization for the 
density in the wind- like case (where p tx r -2 ). As with 
IjMolinari et al.ll2007l) . we assume that n — 1 cm -3 , r\ = 0.2 
and A = 3 x 10 35 cm -1 . 

ipeafc is the time at which the after glow peaks, and 
can b e calculated using Eqn[8] taken from IWillingale et al.l 
l|2007f) : 



peak 



(TrTaY 
V OL a ) 



(8) 



The derived values of t pea k and the initial bulk Lorentz fac- 
tors are shown in Table 2] The calculated values show that 
an ISM type circumburst medium leads to higher Lorentz 
factors, which is consistent w i th th e results shown for a 
larger sample by lEvans et al. | (|201ll . Submitted). By defi- 
nition the afterglow-dominated model peaks at an earlier 
time, which means that the initial bulk Lorentz factor is 
necessarily higher for this CELS6, ELS demonstrated in Table [4] 
It is worth noting, however, that when considering the er- 
rors quoted, the prompt-dominated values cannot be said to 
be distinct from the corresponding values derived from the 
afterglow-dominated fit. 



5 CONCLUSIONS 

In this work, we h ave taken the pulse model of 
iGenet fc Granotl f2 009) and used it in a similar manner 
to IWillingale et al.l ( 20ich to reproduce the prompt light 
curve of GRB 080310. C ombining it with an afterglow model 
l|Willingale et al.ll2007l) . we have tried to produce a simulta- 
neous fit to not only the data from the Swift XRT and BAT 
instruments, but also the available optical and IR datasets. 
The aims behind this work were to establish the origin of 
the observed early-time optical and IR emission, and to at- 
tribute it to either the prompt or afterglow component of 
the GRB. 
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The first conclusion of this work is that a low-energy 
break is required in the spectra of prompt pulses in order to 
fit the optical and near IR flux, regardless of the origin of 
this emission. 

The simplest model considered for the optical and IR 
was to use this low energy break to remove the prompt com- 
ponent entirely from the observed emission. Whilst this suc- 
cessfully recreated the broad scale structure of the optical 
and IR light curves, and the SEDs also appear satisfactory, 
the value of p a to which the fitting tended to was unphyis- 
cal (p a = 5.6) which is inconsistent with that expected for a 
realistic spectrum, such as synchrotron radiation below the 
minimum emitted photon energy (Em) for which p a should 
be ~ or synchrotron self absorption, where p a is expected 
to be | or 2, when the self absorption frequency is above 
or less than E m respectively. However, having looked at the 
one dimensional x 2 surface for p a , we found that the distri- 
bution asymptotes to a better fit at large values of p a . As a 
result we used a more realistic value of p a = §, which didn't 
significantly reduce the quality of the fit. The implications 
of this are that self absorption is a necessary mechanism in 
order to fit the optical and near infrared flux observed, when 
assuming the prompt pulses of the high energy light curves 
do not contribute at early times in the lower energy bands. 

Morphological inconsistencies in the light curves re- 
quired the exploration of an alternative solution. This al- 
ternative was to allow the prompt emission to dominate the 
early times of the optical and IR light curves. An initial 
treatment of the prompt radiation, in which a single value 
for p a was assigned to all pulses, proved insufficient to fit 
the data. Following this, by allowing the pulses to have in- 
dependent values of p a , a fit of similar statistical merit to 
the afterglow-dominated model was obtained (Figures 1101 
and Hip. With this model, an additional break was still re- 
quired for all the prompt pulses, but with a variety of values 
of p a . For five pulses (particularly three of the earliest four) 
steep spectral indices were required, tending to unphysical 
values when fitted. As before, in these instances, we found 
in these instances that it was possible to fix the appropriate 
values of p a to | without significantly altering the quality of 
the fit. This again suggests that self absorption could be an 
important mechanism by which the prompt emission is sup- 
pressed in the optical and infrared regimes. For those pulses 
whose spectrum required a break, but not to the extent of 
pulses 1, 2, 4, 8 or 16, we suggest that the break energy is 
at a value between the optical and X-ray bands, but given 
the degeneracy between p a and E a have not fitted it. These 
pulses could then also have a value of p a = |, but peak at 
an energy nearer to that of the optical bands. 

From the results obtained, it is unclear whether the op- 
tical and IR emission of GRB 080310 originates from cen- 
tral engine or afterglow activity. Neither case accurately de- 
scribes all of the data. The afterglow-dominated model is 
insufficient to describe all of the structure seen in the opti- 
cal light curves, however, the SEDs produced, particularly 
at 350 seconds suggest that the optical and near infrared 
are more faithfully represented by afterglow emission. To 
help discriminate between prompt and afterglow emission 
as a source for early time emission a similar analysis is re- 
quired on a larger sample of GRBs. Bursts in such a sample 
have several important pre-requisites. Firstly, good contin- 
uous optical data are required from very early to late times, 



preferably in several bands, simultaneously. Having such a 
data point from KAIT though, highlights that at the earliest 
times high temporal resolution is required too, as GRBs are 
highly variable during their prompt phases. 

GRBs which will be the best candidates for further anal- 
ysis will be those that contain pulse structure in the optical 
light curves. If these pulses are simultaneous with similar 
structure at higher energies, then it is likely they share a 
common origin. In contrast to this, should the pulses occur 
at markedly different times then an alternative mechanism 
must be found to explain their behaviour. The ideal type 
of burst for this analysis would therefore be one which has 
a long duration as observed in BAT (or Fermi GBM) and 
exhibits strong flaring behaviour after the first hundred sec- 
onds in the X-ray regime. Such times are not only feasible for 
ground based follow-up, but also allow for sufficient tempo- 
ral resolution to discern any features at these lower energies. 
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APPENDIX A: OPTICAL AND NEAR 
INFRARED DATA 

Data from the Swift satellite are publically available. In ad- 
dition to this Tables ET1 - 1 A8l detail the ground-based obser- 
vations taken in the optical and near IR regimes, alongside 
the Sfw/f-UVOT observations used in this analysis. 
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Table Al. All near infrared observations which were calibrated to the i<"-band, showing the central time after trigger of each exposure, 
exposure time (T exp ), filter, instrument, flux and flux errors. The filters quoted are the original filters that observations were taken, 
before conversion to the standard filters used in the later analysis. Fluxes are extinction corrected and any external references are cited. 
A full version of this table is available online. 



Time after trigger (s) 


Texp (s) 


Filter 


Instrument 


Flux OJy) 


Flux error (±/^Jy) External sources 


1650 


47 


K 


PAIRITEL 


1469.20 


101.21 


1723 


47 


K 


PAIRITEL 


1395.98 


95.66 


1796 


47 


K 


PAIRITEL 


1405.27 


103.07 


1868 


47 


K 


PAIRITEL 


1315.23 


96.35 


1941 


47 


K 


PAIRITEL 


1212.27 


92.56 


2013 


47 


K 


PAIRITEL 


1171.55 


92.63 



Table A2. All near infrared observations which were calibrated to the H-band, showing the central time after trigger of each exposure, 
exposure time (T exp ), filter, instrument, flux and flux errors. The filters quoted are the original filters that observations were taken, 
before conversion to the standard filters used in the later analysis. Fluxes are extinction corrected and any external references are cited. 
A full version of this table is available online. 



Time after trigger (s) 


Texp (s) 


Filter 


Instrument 


Flux (/iJy) 


Flux error (±/^Jy) External sources 


185 


36 


H 


REM-REMIR 


820.92 


222.12 


264 


36 


H 


REM-REMIR 


1521.60 


225.43 


354 


36 


H 


REM-REMIR 


1578.70 


250.66 


448 


36 


H 


REM-REMIR 


942.54 


201.13 


598 


86 


H 


REM-REMIR 


1122.79 


131.21 


1650 


47 


H 


PAIRITEL 


1174.00 


65.33 



Table A3. All near infrared observations which were calibrated to the J-band, showing the central time after trigger of each exposure, 
exposure time (T exp ), filter, instrument, flux and flux errors. The filters quoted are the original filters that observations were taken, 
before conversion to the standard filters used in the later analysis. Fluxes are extinction corrected and any external references are cited. 
A full version of this table is available online. 



Time after trigger (s) 


Texp (s) 


Filter 


Instrument 


Flux GuJy) 


Flux error (±^tjy) External sources 


1650 


47 


J 


PAIRITEL 


961.77 


44.42 


1723 


47 


J 


PAIRITEL 


903.55 


41.16 


1796 


47 


J 


PAIRITEL 


932.98 


47.86 


1868 


47 


J 


PAIRITEL 


879.73 


43.09 


1941 


47 


J 


PAIRITEL 


928.43 


44.19 


2013 


47 


J 


PAIRITEL 


887.46 


43.31 



Table A4. All optical observations which were calibrated to the /-band, showing the central time after trigger of each exposure, 
exposure time (T exp ), filter, instrument, flux and flux errors. The filters quoted are the original filters that observations were taken, 
before conversion to the standard filters used in the later analysis. Fluxes are extinction corrected and any external references are cited. 
A full version of this table is available online. 



Time after trigger (s) T exp (s) Filter Instrument Flux (/^Jy) Flux error (±/iJy) External sources 



384 60 V P60 650.01 69.31 Cenko et al. (2009) 

642 60 V P60 571.37 60.92 Cenko et al. (2009) 

899 60 V P60 587.38 33.37 Cenko et al. (2009) 

1246 120 i' P60 632.30 35.93 Cenko et al. (2009) 

1448 60 i' P60 620.76 23.30 Cenko et al. (2009) 

1541 45 I c SMARTS-ANDICAM 676.00 32.00 Kann et al. (2010) 
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Table A5. All optical observations which were calibrated to the ij-band, showing the central time after trigger of each exposure, 
exposure time (T exp ), filter, instrument, flux and flux errors. The filters quoted are the original filters that observations were taken, 
before conversion to the standard filters used in the later analysis. Fluxes are extinction corrected and any external references are cited. 
A full version of this table is available online. 



Time after trigger (s) 




Filter 


Instrument 


Flux (ply) 


Flux error (±/^Jy) 


External sources 


57 


30 


unfiltered 


KAIT 


29.70 


9.45 




162 


20 


unfiltered 


KAIT 


571.23 


10.62 




171 


20 


R 


Super-LOTIS 


218.11 


39.33 


Milne & Williams (2008) 


195 


10 


R 


Super-LOTIS 


634.87 


48.55 


Milne & Williams (2008) 


224 


36 


R 


REM-ROSS 


683.41 


161.25 




254 


20 


unfiltered 


KAIT 


781.29 


7.23 





Table A6. All optical observations which were calibrated to the V-band, showing the central time after trigger of each exposure, 
exposure time (T exp ), filter, instrument, flux and flux errors. The filters quoted are the original filters that observations were taken, 
before conversion to the standard filters used in the later analysis. Fluxes are extinction corrected and any external references are cited. 
A full version of this table is available online. 



Time after trigger (s) 




Filter 


Instrument 


Flux OJy) 


Flux error (±/^Jy) External sources 


84 


10 


V 


UVOT 


246.60 


140.41 


113 


10 


white 


UVOT 


115.97 


43.37 


123 


10 


white 


UVOT 


62.36 


40.39 


133 


10 


white 


UVOT 


102.54 


43.27 


143 


10 


white 


UVOT 


215.19 


50.65 


153 


10 


white 


UVOT 


321.34 


57.17 



Table A7. All optical observations which were calibrated to the 6-band, showing the central time after trigger of each exposure, 
exposure time (T exp ), filter, instrument, flux and flux errors. The filters quoted are the original filters that observations were taken, 
before conversion to the standard filters used in the later analysis. Fluxes are extinction corrected and any external references are cited. 
A full version of this table is available online. 



Time after trigger (s) 


Texp (s) 


Filter 


Instrument 


Flux (fjjy) 


Flux error (±^tjy) External sources 


690 


10 


b 


UVOT 


451.91 


87.74 


843 


10 


b 


UVOT 


382.97 


82.63 


1465 


20 


b 


UVOT 


406.60 


79.25 


5247 


199 


b 


UVOT 


213.69 


15.71 


6682 


199 


b 


UVOT 


139.11 


14.72 


23481 


906 


b 


UVOT 


40.52 


4.67 



Table A8. All optical observations which were calibrated to the u-band, showing the central time after trigger of each exposure, 
exposure time (T exp ), filter, instrument, flux and flux errors. The filters quoted are the original filters that observations were taken, 
before conversion to the standard filters used in the later analysis. Fluxes are extinction corrected and any external references are cited. 
A full version of this table is available online. 



Time after trigger (s) 


T eX p (s) 


Filter 


Instrument 


Flux (fjjy) 


Flux error (±/xJy) External sources 


670 


20 


a 


UVOT 


143.74 


35.38 


823 


20 


u 


UVOT 


114.74 


33.43 


1440 


20 


u 


UVOT 


261.99 


48.14 


5041 


200 


u 


UVOT 


152.79 


11.27 


6477 


200 


U 


UVOT 


102.75 


10.28 


22568 


906 


U 


UVOT 


23.23 


2.85 



